Maximally entangled gapped ground state of lattice fermions 
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Entanglement between the constituents of a quantum system is an essential resource in the im- 
plementation of many quantum processes and algorithms. Indeed, universal quantum computation 
is possible by measuring individual qubits comprising highly entangled 'cluster states.' In this work 
it is shown that the unique gapped ground state of non-interacting fermions hopping on a specially 
prepared lattice is equivalent to a cluster state, where the entanglement between qubits results 
solely by fermionic indistinguishability and antisymmetry. A deterministic strategy for universal 
measurement-based quantum computation with this resource is described. Because most matter 
is composed of fermions, these results suggest that resources for quantum information processing 
might be generic in Nature. 



I. INTRODUCTION 

In measurement-based quantum computation 
(MBQC) [1], an algorithm proceeds entirely by making 
projective measurements of successive particles (qubits 
or qudits) comprising some highly entangled 'resource 
state.' While two-dimensional cluster states are known 
to be universal resources for MBQC |2j, it has been 
proven that they cannot be the unique ground states 
of any qubit Hamiltonian involving at most two-body 
interactions @. Because two-body interactions are 
generic in Nature, this limitation has driven a search 
for other universal resource states that might exist as 
the ground state of a physical system. An alternative 
representation of MBQC using Matrix Product States 
(MPS) 4] in one dimension and its extension to higher 
dimensions using Projected Entangled Pair States 
(PEPS) indicates that there are many states with 

much less entanglement than the cluster states that 
are nevertheless universal resources 043- Unfortu- 
nately, the fraction of useful resource states is known 
to be exponentially small in the number of physical 
qubits [Hi El- Worse, it is presently unknown if any of 
these states is the unique, gapped ground state of any 
physically realizable Hamiltonian. By relaxing various 
assumptions, it is possible to define accessible Hamil- 
tonians whose ground states are universal resources for 
MBQC. One example is a system of spin-5/2 qudits on 
the hexagonal lattice with two-body interactions, where 
the ground state is unique and gappe d 11211 . Another is 
the AKLT model of an antiferromagnet[13| with spin-3/2 
qudits on the hexagonal lattice, which is gapped but 
degenerate flil. [l5j. 

An unrelated strategy for MBQC was devised for in- 
distinguishable free fermions using beam splitters and 
charge detectors [l6| in analogy with non-deterministic 
linear optical quantum computation [l7j . In fact, clus- 
ter states can be dynamically generated using this ap- 
proach (l8j . The relationship between free fermions and 
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entanglement has subsequently been the subject of ex- 
tensive investigations [19|, |20j. More recently, the re- 
lationships between fermions and quantum information 
are being explored pll - t^j ] . In particular, it has been 
shown that quantum circuits composed only of match- 
gates [IE Hl| (which are closely related to non- interacting 
fermions) can generate entanglement yet are efficiently 
classically simulatable unless supplemented by an appar- 
ently trivial resource, the unentangling SWAP gate |22| . 

In this work, the ideas above are combined to show that 
the non-degenerate and gapped ground state of a non- 
interacting fermionic lattice Hamiltonian is equivalent to 
a cluster state comprised of qubits, under the conditions 
that any site has at most one fermion and the lattice is 
exactly half-filled. The qubit registers then correspond to 
the occupation of lattice sites. The entanglement is a di- 
rect consequence of fermionic indistinguishability and an- 
tisymmetry. Measurements correspond to detecting the 
presence or absence of a fermion in a lattice site, equiva- 
lent to a restriction to the computational basis. In order 
to simulate a universal quantum computer, one must ro- 
tate the measurement from the computational Z basis 
into the X or Y basis via the application of a Hadamard 
gate if or a \[X operation (here X, Y, and Z repre- 
sent the Pauli matrices). Surprisingly, in the fermionic 
model these gates cannot be directly implemented be- 
cause of the intrinsic entanglement; rather, an additional 
lattice must be introduced that allows two fermions en- 
coding different qubits to interact, thereby cancelling the 
effects of the fermionic antisymmetry. Fermionic systems 
therefore invert the usual quantum computing paradigm, 
in that entanglement is intrinsic but general single-qubit 
operations require particle interactions. 

The manuscript is organized as follows. A review of 
cluster states and the measurement-based model for uni- 
versal quantum computation using them is briefly re- 
viewed in Sec. [TTJ Sec. IIIII introduces the main fermionic 
model, and demonstrates that in principle the ground 
state is a universal resource for MBQC. The implementa- 
tion of a universal set of quantum gates by measurements 
is discussed in Sec. lIVl where it is shown that interactions 
and additional potentials are required. The implications 
of the results are summarized in Sec. [Vj 
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II. REVIEW OF CLUSTER-STATE QUANTUM 
COMPUTATION 



The stabilizers of a given state are the operators 
Si for which one satisfies Si\ip) = for each i. If 
is a stabilizer state, then its Hilbert space dimension ex- 
actly coincides with the number of stabilizers Si. A useful 
class of stabilizer states are known as graph states, where 
qubits are represented by vertices and maximal two-qubit 
entanglement between two qubits is represented by an 
edge connecting the vertices. For a generic graph state 
\g) on N qubits, the N generators for the 2 Ar -dimensional 
group of stabilizers are given by the separable operators 
Si — Xi 0^^) Zj, where Af(i) signifies the neighbour- 
hood of site (qubit) i, and where Xi, Yi, and Zi denote 
Pauli matrices acting on qubit i. A Hamiltonian that is 
guaranteed to have \g) as the lowest energy state (with 
eigenvalue —Nt) is therefore 



\— g ) = -J= (l — e l9 ); alternatively, the operation 

HRz(0) is first applied to the qubit, and then it is 
measured in the computational basis. In both cases, 
the result of measuring qubit 1 of a linear cluster state 
with N qubits is |m)i ® [X m HR z (9)} 2 \+ 9 V f - 1) )n,3,...N, 
where m = {0, 1} is the measurement outcome and the 
subscript is the qubit index. The presence of the X 
gate reflects the probabilistic nature of the measurement 
outcome. Three successive measurements with three 
independent angles (adapted to previous measurements 
because of the possible presence of an X gate) suffices 
to telcport an arbitrary single-qubit gate in ID cluster 
states; a universal set of gates results in 2D cluster 
states 



III. ENTANGLEMENT IN LATTICE FERMIONS 



N 



H 



(i) 



For cluster states, in which the qubits are arranged on 
a regular lattice, the coordination number is two and 
four for one-dimensional (ID) and two-dimensional (2D) 
lattices, respectively. This implies that the Hamiltonian 
above consists of three-spin and five-spin operators in 
ID and 2D, respectively. In fact, this number has been 
proven to be a minimum, i.e. there is no Hamiltonian 
having all terms with fewer spin operations that can have 
\g) as a (non-degenerate) ground state 0]. 
In the one-dimensional case, the cluster ground state 

is 



N-l 



\C 



N 



n rz << 



V 



(2) 



where CZij = diag (1,1,1,-1) is the maximally entan- 
gling controlled-phase gate and |+) = (|0) + |l))/\/2 is 
the +1 eigenstate of the X operator. In this representa- 
tion, it is evident that the cluster state is highly entan- 
gled: each qubit is maximally entangled with its neigh- 
bor. The excited states correspond to Ylf = i Z^ J \Cn) , 
where kj <E {0, 1}, and have energies (— N + £)t with 
I € {1, . . . , 2N}. The energy gap to the first excited state 
is therefore r, constant in the thermodynamic limit. 

To perform universal quantum computation, it suffices 
to measure qubits in successive columns of the 2D lattice 
in the HRz{9) basis, where H and Rz{9) = exp (— i%Z\ 
are the Hadamard operator and rotation by an angle 9 
around the Z-axis, respectively: 



-iB/2 

e w / 2 



This corresponds to measurements in a basis 
spanned by the vectors |+e) = 4j (l + e~ ie ) and 



A. Fermionic model 

The fermionic model employed in this work is inspired 
by experiments with ultracold atomic gases. Numerous 
groups have confined both bosonic and fermionic atoms 
in periodic potentials formed from interfering lasers, 
known as optical lattices [13, [Hf . These lattices can con- 
fine different spin components of the same bosonic species 
in independent, overlapping optical lattices 29,[30j]. The 
structure of these lattices is highly flexible [31( , and pe- 



riodic double-well potentials for bosons have been been 
the subject of particularly intense investigations [32h36| . 
In another recent development, an ultracold gas consist- 
ing of a mixture of three different spin components of 
the fermionic atom 6 Li has been realized in an optical 
trap [37], l38j . Ultracold fermionic gases of both 40 K 
and 6 Li have been loaded into periodic optical lattices 
with highly controllable geometries and trapping parame- 
ters [39|-[42[ . It is reasonable to assume that these various 
techniques could be combined in the future, with multi- 
component fermions in three overlapping double- well op- 
tical lattices, each confining a different spin component. 

Consider the potentials depicted in Fig. 1(a), corre- 
sponding to three overlapping one-dimensional two-site 
lattices. It is assumed that each sublattice can confine 
only one particular spin projection of a fermion with 
pseudospin s — 1 (corresponding to real spin s > 3/2 
with all but three spin projections frozen out); the spin 
projections will be labeled f> °, an d i- The potential 
shown as a black curve is assumed to confine only the 
spin o, and could be obtained using two potentials with 
amplitudes V\. and V^o and with one wavelength A half 
as large as another: 



V (x) = ^i, o C0S 2 ( — ) + V2, o C0S 2 



A 



Likewise, the potentials seen by the f and J. spin pro- 
jections (green and red curves in Fig. HJa), respectively) 
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FIG. 1. (Color online) Depiction of the fermionic leapfrog 
model. Each spin projection experiences a unique spatially- 
dependent potential energy, shown as the black, red (dark 
gray), and green (light gray) curves in (a). Assuming that 
the higher potential barriers are effectively impenetrable, one 
obtains the hopping model in (b), with three overlapping 
double-well sublattices. A two-dimensional extension of the 
one-dimensional model is shown in (c). 



overlap in practice. For example, a o-projection fermion 
in site 1 in Fig. 1(b) needs to be antisymmetrized only 
with fermions in sites 2, 3, and 4; the tunneling barrier 
between sites 4 and 5 is assumed to be so large as to 
prevent any tunneling for o-spin projection atoms. 

The relative strengths of the tunneling amplitudes are 
not important for the ground state, as long as they are 
all non-zero. One may assume that the system is given 
sufficiently long to equilibrate compared to the tunneling 
rates. Thus one may choose t q = t-j- = rj, = r with- 
out loss of generality. The Hamiltonian ^ can be fur- 
ther simplified by noting that fermions with a given spin 
projection are constrained to their respective sublattices. 
The spin-projection index can therefore be suppressed, 
corresponding to spinless (fully polarized) fermions. One 
then obtains the equivalent 'leapfrog' Hamiltonian 

h = -tY, (flj+ihj+i + fh+ihj+i) , (4) 

which is the central model of the current work, and is 
depicted graphically in Fig. 1(b). Though the Hamilto- 
nian appears to involve third-nearest neighbor hopping, 
it is in fact obtained by nearest-neighbor hopping on an 
array of three overlapping two-site sublattices. 



would be generated by wavelengths that were twice as 
long as for projection o, and offset relative to each other: 



V\(x) = Vi,t cos 2 



ir(x — A) 
A 



V i (x) = V 1 , l cos 2 ( K —)+V 2 , i < 



V2,f COS 

2 f nx \ 
)S 

V2A/ 



ir(x — A) 
2A 



The relative height of the barriers is then set by the ra- 
tio V\, a jV2.G- Within the tight-binding approximation, 
assuming the atoms can only tunnel through the lower 
barrier of a given sublattice, the Hamiltonian becomes 



N-l 



H = - ( T t/8i+2,t/8i- 1 -t + T ofl ]+4: ,ohj+l,o 



3=0 



+n/si+64/8i+84 + H - c -) 



(3) 



where fj (fi,a) creates (annihilates) a fermion with spin 
projection a at site i. The r CT correspond to hopping 
amplitudes, the exact values of which are unimportant 
for the present work. 

There is no tunneling permitted between different sub- 
lattices, but nevertheless fermions with different spin 
projections must anticommute if the sublattices over- 
lap: if a fermion with one spin projection tunnels past a 
fermion with another spin projection, the total wavefunc- 
tion picks up a negative sign. This is because it is im- 
possible to distinguish if the fermions interchanged their 
spin labels (switching places) during the transit, or re- 
mained the same. At the same time, it is unnecessary to 
impose antisymmetrization on fermions that can never 



B. Cluster ground state 

Consider only the case where each spinless sublattice 
is exactly half-filled: each connected double-well indexed 
by j contains exactly one fermion, ri2j+i + n2j+4 = 1- 
Mathematically, this could be accomplished by applying 
a suitable projector to the Hamiltonian (U). In an experi- 
ment with ultracold atoms, such a configuration could be 
obtained by detecting and selectively moving or removing 
single atoms [43T - |47j . 

Making use of the Jordan- Wigner transformation 



1 



fe=r 



/ 



1 



3-1 



\Z k ®(X i -iY j ), 



(5) 



fe=r 



the Hamiltonian (j3J may be expressed in the spin repre- 
sentation: 



7V-1 



H — —— ^ Z 2 j+2Z2j+3 (X2j+lX2j+4 + Y 2 j + lY2 



2j+4) , 



3=0 



which is four-local in marked contrast to that in the 
fermion representation. Because of the half-filling restric- 
tion, the ground state only involves states with Hamming 
weight TV and therefore only involves combinations of the 
encoded registers 

|Q J .) = |l 2i+1 2j - +4 ) = 4. +1 |O); 

|l J ) = |0 2 , +1 l 2j+4 )=4 +4 |O), 
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where the state \0) represents the particle vacuum. In 
other words, the \0j) register has the single fermion on 
the first site of a given double- well indexed by j, while 
the |1_-) register corresponds to the fermion in the second 
site. 

One can introduce an encoded basis where the op- 
erators become (X 2 j+iX 2 j+4 l + Y 2 j+iY 2 j+i) /2 = X_p 
hj+iZ 2 j + 4 = Z_j, and Z 2 j + iI 2 j + 4 = -Z_ rj (recall 
the half-filling condition), i.e. Xj\Qj) = \lj) and vice 
versa, Z^Oj) = |0^), and Z^lj) = Using 

%2j+2 = i2(i-l)+1^2(j-l)+4 = an d %2j+3 = 

■^2(j+i)+i^2(j+i)+4 = the Hamiltonian ^ in the 

encoded basis becomes 

JV-l 

H = t ^2 Kj-i2LjZ j+1 - tZ n _ 1 X n , (7) 

J=l 

with Z_ = Applying the unitary transformation 

Z_j for all j e { 1 , AT — 1 } performs the maps 2Lj 
—X_j, transforming ([7]) into exactly the same Hamilto- 
nian as the one-dimensional cluster-state Hamiltonian on 
N qubits, Eq. ([T]). The ground state of Hamiltonian ([7]) 
is therefore immediately obtained by comparing with the 
usual cluster state, Eq. ([2]): 

N-l 

ig. s .)= hz 3 cz^ +1 \+® n ), 

where |+) = (|0) + |1)) /v2, and the encoded entangling 
gate is defined as QZ_ j j + i = diag (1, 1, 1, —1) in the space 
spanned by the vectors \0jQj +1 ), \Q.jl j+ i), |ljQ,,- +1 ), and 
lljlj+i)- It is clear that the ground state is equivalent 
to the one-dimensional cluster state, where each qubit is 
maximally entangled to its neighbor. This entanglement 
is a direct consequence of fermionic antisymmetry and 
indistinguishability. The excitations correspond to Z^ 
operations on the ground state, which cost energy r. The 
energy gap is therefore r, independent of system size. 

C. Example: two qubits 

Consider the consequences of measuring the position of 
one of the two fermions in a two-qubit fermionic cluster 
state. Suppose that the entire system consisted of sites 1, 
3, 4, and 6 in Fig. 1(b). In the absence of inert fermions 
in the unpaired sites 2 and 5, these sites can be ignored, 
and the ground state of the Hamiltonian (|4]) is maximally 
entangled: 

|g-B->-^(/l f + /|) (fl + ft) P) 1,3,4,6 

= \ (|1100) + |1001) - |0110) + |ooii)) 1)3i4)8 

= Zi^i 2 l++)- (8) 

In the Fermi-occupation notation, the ground state ap- 
pears to be separable; however, this ignores the fact that 



fermions on sites 3 and 4 must anticommute if there is 
any chance that these fermions can interchange. Con- 
sequently, defining the fermionic occupation states in 
the natural ordering by increasing label |01 10) 1,3^^ = 
/j/llO) 1,3,4,6 1 a negative sign appears on the second line 
above. If the fermions on site-pairs {1,4} and {3, 6} could 
never overlap (i.e. there was an impenetrable barrier be- 
tween them) , then the state would be truly separable and 
no explicit antisymmctrization would be required. 

Measurement of the first qubit in encoded state jOj) 
corresponds to the application of the projector = 
|1i0 4 ><1i04| = f\\0){0\h, yielding 

(Mi) 1 |g.s.) = -^[|1100) + 11001)] = |0+). 

Likewise, applying the other projector (^2)1 = 
|0il 4 )(0il 4 | = fl\0)(0\U yields the state 

(Ma)! |g.s.) = i [-10110) + 10011)] 

= -£i|l+) = -|l-). 

The dependence of the second encoded qubit's output 
state on the outcome of the first encoded qubit's mea- 
surement is a direct consequence of the entanglement of 
the original fermionic state. 

It is important to note, and easy to verify, that this 
result doesn't depend on how one labels the sites. If the 
two double-sites were instead labeled {1,3} and {4, 6}, so 
that no negative sign appeared in the ground state ((5), 
the same states would still be obtained after measure- 
ment. This is because the result of the measurement op- 
eration on the double-site {1,3} explicitly depends on 
the possible existence of a fermion at the intervening 
site 4. If instead the two double-wells were disconnected 
(non-overlapping), then there would be no intervening 
site with which to anticommute, and therefore no de- 
pendence of the output on the measurement outcome. 
In short, the entanglement arises from the fundamen- 
tal requirement that fermions on the sites 3 and 4 must 
anticommute when there is the potential for these wave- 
functions to overlap. 



IV. UNIVERSAL MEASUREMENT-BASED 
QUANTUM COMPUTATION WITH LATTICE 
FERMIONS 

A. Single-qubit operations with fermions 

In order to perform universal gate teleportation using 
the encoded cluster state, Eq. (10) in the manuscript, it 
is important to derive expressions for the encoded X_, Y_, 
and Z-rotations, Hadamard operator _ff, and measure- 
ment operators that one would need to apply in order to 
teleport gates with fermions. The encoded Z-rotation is 
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straightforward: 



be 



Rz{0)j 



exp 



cos 



-iZ A 



exp 



-iZ- 



— i sin 



(1 - 2n 2j+4 ) , 



where 1 — 2ri2j+4 = Z_j is the encoded Z operator. Be- 
cause U2j+A — 1 — ?^2j+i an equivalent expression would 

I 



cos 



i sm 



(1 - 2n 2 j+i) 



Thus, i? z ((9) = e- 10 ' 2 if the fermio n is in the first site 
and Rz{0) = e l6 < 2 if it is in the second, as expected from 
the qubit encoding. To obtain the encoded X-rotation, 
one can make use of the relations 



X j = (1 - 2n 2i+2 ) (1 - 2n 2j+a ) 
x {jtj+ihj+i + ftj+ifaj+i) 



(9) 



and 2Lj = n 2]+ i (I - n 2 j+4) + njy+4 (1 - n 2 j+i) = 
U2j+x + ^2j+4 = 1 when each two-site lattice has exactly 
one fcrmion. One then finds 



Rx (<% = ex P (-iXj^j = cos ^0 - isin f (1 - 2n 2j+2 ) (1 - 2n 2j+3 ) (jv+Jv+i + flj+ihj+i) 



One can write Y_ rj = —iZ_jX_j = iZ 2 j + iXj, so that the encoded y-rotation is 



R Y {6)j = exp (-iYj^j = cos f^j + sin O^J (1 - 2n 2j -+i) (1 - 2n 2 j +2 ) (1 - 2n 2 j +3 ) [flj+ihj+i + flj+ifci+t) 



Finally, the encoded Hadamard gate is H_j = Z_jR Y (—ir/2)j = — (1 — 2n 2 j_|_i) Ry(— n/tyj'- 

1 



-1 + 2n 2j -+i + (1 - 2n 2i+2 ) (1 - 2n 2j+3 ) \fl j+ ihi+i + fl j+ j2j+4 



(10) 



Note that in marked contrast with the qubit represen- 
tation, in the fermion representation the 2Cp Y_ji an d 
H^j operators acting on encoded qubit j (involving the 
two sites labeled by 2j + 1 and 2j + 4) are intrinsically 
non-local, since they depend explicitly on the fermion 
density at the intervening sites 2j + 2 and 2j + 3. Nev- 
ertheless they are fully separable (i.e. not entangling), 
as will be shown explicitly in the example below. One 
might ask why they are nevertheless non-local. Suppose 
that the (1 — 2n 2 j+ 2 ) (1 — 2n 2 j+z) term were dropped in 
expression ^ for the encoded X_- operator. The op- 
eration would then naively correspond to an unencoded 
Xj operation, but only in the absence of fermions be- 
tween sites 2j + 1 and 2j + 4. If there were in fact one 
fcrmion occupying either of the intervening sites 2j + 2 
or 2j + 3, then hopping between 2j + 1 and 2j + 4 would 
introduce undesired negative signs to certain matrix el- 
ements. Of course, these are exactly the factors that 
yield the ground-state entanglement in the first place! 
But they are unwanted in the implementation of X, Y, 
and H operations, and so need to be cancelled in the 
encoded versions X, X_, and H_ if necessary. Consider 
two encoded qubits, for example. The basis is spanned 



by the vectors |1100) = flf 3 \0 4 ), |1001) = ///eTO, 
10110) = ftflpi), and 10011) = flftpi), the encoded 
Hadamard gate (fTO)) acting on the first encoded qubit is 




V2 



(X + Z) 



(11) 



which is manifestly separable in matrix form even though 
it appears not to be in the fermion representation (|10[) . 

The measurements are made in the local fermion basis, 
so that 



(Mi),- 



|l2j + l0 2 j+4) (l2j + l0 2 j-(-4| 



"-2j + lU 2 j+4/ 
= f 2 j + lp2)2j+l,2 3 +4(0 2 \f 2j+ l 



and 



{M 2 ) 3 = |1),(1| = |0 2j+1 1 2j+4 )(0 2j+1 1 2j+4 | 



4 +4 P 



2i+ 4l^2/2i+l,2i+4 



(02 1/: 



2j+A- 



In other words, the qubit projectors correspond simply 
to observing the site indexed by 2j + 1 or 2j + 4 to de- 
termine if a fermion is present; if so, it is projected into 
the encoded state |0) or |1_), respectively. 



B. Example: two-qubit gate teleportation 

It is instructive to consider how one would implement 
one gate teleportation on a two-qubit fermionic cluster 
state. Suppose that the entire system consisted of sites 
1, 3, 4, and 6 in Fig. 1(b) as in the two-qubit exam- 
ple in the previous section, with ground state given by 
Eq. ©. Suppose that all sources of energy exchange be- 
tween the system and the environment are eliminated, 
other than those associated with applied rotations and 
measurements, so that \ip) = |g.s.). In order to teleport 
an arbitrary single-qubit gate in a ordinary ID cluster 
state, one needs to perform a series of measurements in 
the HRz(O) basis. The first qubit is to be rotated first 
by the encoded gate Rz(@)i, then by H_ x , and then mea- 
sured in the computational basis spanned by the states 
|0> = fl\0) 1A and |1) = fl\0) 1A . Evidently the applied 
rotations on the first qubit will push \ip) away from the 
ground state, though the fermions in all the other double 
wells remain in eigenstates of the original Hamiltonian. 
As a consequence, the state of the first qubit will oscil- 
late at a frequency th/H. All states and operations can 
therefore be considered only at 2tt /t\a intervals. This as- 
sumption is not essential, however, as discussed below in 
Sec. UVDl 

A j^-rotation on the first encoded qubit yields 

W) = EzVM = \ (/i + e W fl) (/] + A) I<?>1,3,4, 6 
= \ [|1100) + |1001) + e* e (-10110) + |0011))] 1)M)6 , 

ignoring an overall phase e _ie / 2 ; subsequent application 
of an encoded Hadamard operation gives 

W) = K,W) = -±= [(1 - e ie ) (|1100) + |0011)) 
+ (1 + e ld ) (|1001) + |0110))] li3A6 . 

Measurement of the first qubit in encoded state |0 X ) cor- 
responds to the application of the projector = 
|li0 4 )(li0 4 H/ 1 t |O)(O|A: 

{Mx) x W) = \ [(1 - e ie ) |1100) + (1 + e* e ) |1001)] 

= ||0x> [(l-e t9 )|0 2 ) + (l + e ^) |1 2 >] 
= K 2 Z 2 R Z (9) 2 \0+). 

Likewise, applying the other projector (A4 2 )i — 
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FIG. 2. (Color online) A subgraph of the two-dimensional 
leapfrog model, with the two overlapping sublattices corre- 
sponding to horizontal red (gray) and black links. 

|0il 4 )(0il 4 | = fl\0){0\f 4 yields the state 

(M 2 \ |V") = \ [(1 + e* e ) |0110) + (1 - e* e ) |0011)] 

= liii) [(l + e ie )|Q 2 ) + (W 9 ) |1 2 )] 

= 2L 2 H 2 Z 2 R z (9) 2 \l+). 

Thus, the gate teleportation using fermions proceeds ex- 
actly as in the usual spin-qubit case, though with an ad- 
ditional overall Z_ operator reflecting the original ground 
state ©. 

C. Two-dimensional cluster state 

Though the Jordan- Wigner transformation is not 
uniquely defined in two dimensions, the structure of 
the qubit encoding carries naturally over to dimensions 
higher than one. A two-dimensional version of the 
leapfrog Hamiltonian is shown in Fig. [He). Two in- 
terpenetrating two-site lattices (shown in grey) that are 
staggered relative to each other are superimposed on an 
array of one-dimensional leapfrog chains, similar to those 
discussed above. Each double-well of the new lattices 
serves to connect adjacent one-dimensional chains, so the 
sites are situated between a given connected pair in the 
chains. Again, only a single fermion is allowed to occupy 
any given pair of sites of the new lattices. The labelling 
of the sites can be chosen in any convenient way. When 
writing out the full ground state in terms of the fermion 
field operators, one need only keep in mind which sites 
interleave others. 

In order to demonstrate that the ground state of the 
fermionic leapfrog Hamiltonian is computationally uni- 
versal, it remains to be shown that the presence of 
the vertical links has the effect of entangling the two 
horizontal chains. Consider for this purpose a two- 
dimensional subgraph of the full lattice configuration 
shown in Fig. QJc), in which two encoded two-qubit clus- 
ter states are coupled together by a single vertical link 
as shown in Fig. [5] The sites on the top row are labeled 
1 through 5 with site 2 the upper partner of the verti- 
cal double well, and the bottom three sites are labeled 5 
through 10 with site 7 the lower partner of site 2. The 
Hamiltonian corresponding to this subgraph is 

ff 2D = -r (/i/i + Ah + Ah + Ah + Aoh + h.c.) . 
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It is simple to verify that the ground state is 

|g-) 2D = -^ f (/?+/!) (fl + fl) (ti + ff) 

[fl + fS) (.4 + flo) P). 



application of the additional two-site hopping Hamilto- 
nians 



H a (ri4, V u Vi) = m [ flfi + flh + Vim + V A n A 



If the vertical link (sites 2 and 7) truly entangles 
the two horizontal chains, then measuring the fermions 
in all the double wells labeled by {1,4}, {6,9}, and 
{2, 7} should yield entanglement between the two en- 
coded qubits associated with the double wells labeled by 
{3,5} and {8,10}, even though there is no explicit link 
remaining between them. The {1,4} and {6,9} sites are 
measured in the encoded X basis, corresponding to the 
application of an encoded Hadamard gate followed by a 
computational-basis measurement. The {2, 7} sites are 
measured in the encoded Y basis, i.e. the {2, 7} sites are 
first transformed by \[X and then measured in the com- 
putational basis. The result doesn't depend on the order 
in which these operations and measurements are made, 
because only Pauli byproduct operators are teleported. 
Ignoring overall phases the output is 

|out) = |mim2m 3 )i,2,4,6,7,9 

<8 Uo-ffi \fZ~ x Ho yfzlCZ^ 1 ++)! 

L J 3,5,8,10 

where nii — {0, 1} are the measurement outcomes and 

are local encoded Pauli byproduct operators correspond- 
ing to the eight possible outcomes. Thus, the result of 
these measurements is to maximally entangle the two 
qubits |+)3,5 and |+)s,io- If the two qubits encoded on 
sites {1,4} and {6,9} were originally in the states \if) ) 
and \ip 2 ), then the result would be the same as that above 
except with |++) 3 ,5,s,io -> l^^) 3,5,8,10- 

D. Practical Considerations 

The two-qubit example demonstrates how one- 
dimensional fermionic cluster states can be used in order 
to perform universal gate teleportation. Nevertheless, it 
is not apparent how the various encoded operations — 
the R z (9), H_, R x (tt/2) = \fX. and measurements - 
would be implemented in practice. Of particular con- 
cern is that the encoded Hadamard gate and the X_ and 
Y_ rotations are explicitly non-local in the fermion rep- 
resentation, as shown above [cf. Eqs.(j9|) and (fT0|)]: all 
require local hopping on a given pair of sites in a sub- 
lattice with an amplitude dependent on the occupancy 
of intervening qubits on different sublattices. As shown 
below, it is possible to implement all encoded operations 
but only if one allows the fermions to explicitly interact. 

Consider again the two-qubit example discussed above. 
All local gates on the first qubit can be performed by the 



and 



(12) 



H 



b (T34,5, V 3 , Vi) = T34 \jl >rT h,o + fl, a U,a + 5™3,<7™3,a' 

+V 3 n 3 ^ + Vmi, a ) , (13) 



for particular choices of parameters and times. The 
Hamiltonian (fT2|) is the same as Eq. ((4]), with the ad- 
dition of potential energies V\ and V4 on the two sites of 
the first double-well. The spin indices a and a' in Hamil- 
tonian p3[) represent the fermionic spin projection for 
encoded qubits 1 and 2, respectively. Thus, the fermion 
encoding qubit 1 is allowed to tunnel between sites 4 and 

3 and possibly interact with the inert fermion on site 3. 
A fermion may or may not be present on site 3 (if not, 

then it occupies site 6), so the Hamiltonians (jl"2j) and 
(H2J) act on two independent two-dimensional blocks. For 

H a , the states in the first block are j/J/^O) , flft\0}\, 

while those in the second block are j/j \0) ,f\\0)\. In 

the first block, a fermion hopping between sites 1 and 

4 must cross the fermion in site 3, acquiring a negative 
sign in the process; the sign change is absent in the second 
block. But all single-qubit operations on the first qubit 
must be truly local: the resulting gate ex-p(-itH) cannot 
depend on the state of the second qubit, i.e. the possible 
presence of a fermion on site 3. For H , the two states in 

the first block correspond to [f^fljO) , f{ a , f{ a P)) , 

while those in the second block are |/3 CT |C) , /lo-IC)}- 

No fermions cross in either block, but additional phases 
can result from the interaction in site 3 of a fermion with 
spin a and one with spin a' . 

The most straightforward gate is an arbitrary Z- 
axis rotation on the first qubit, which corresponds 
to turning on H a for a given time: H z {— V A i)\ — 
cxp \—itH a (0, 0, V4)] for either block, ignoring an overall 
unimportant phase. In practise, a strong repulsive poten- 
tial would be applied to the barrier between sites 1 and 
4 to turn off tunneling, and another potential applied to 
site 4. Next consider the implementation of the encoded 
Hadamard gate H_ 1 . Applying H a (tu, — 1 + v2, 1 + v2j 
for a time t — 7r/v^Ti4 yields H_ 1 in the first block 
and Z_iM_ l Z_ l in the second. This gate is equivalent 
to C 'Z 12 H \ C 'Z 12 which is not the desired Hadamard 
gate. Likewise, the application of H a (ti4,0,0) for a 
time t = 7t/4ti4 yields C Z 12 R x {-k 1 2)y C Z l2 rather than 
R x (n/2)i. In order to obtain the desired gates, one can 
apply Hf, both before and after the H a above. Choosing 

Hb (t 34 , -^=, -^j for a time t = V3Tr/2T 3i yields Jj 

in the first block and Z_ x in the second, i.e. the gate CZ 12 . 
All parameters in H a are assumed to be zero during the 
operation of H . Because Hi, also applies a possible phase 
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to the first register of the second qubit, a repulsive a' po- 
tential should be applied between sites 4 and 6 during its 
application and removed thereafter. The majority of the 
state thereby remains in the ground state of the primary 
governing Hamiltonian. 

It is interesting to note that any attempt to perform 
a local operation U_ 1 within the inital non-interacting 
fermionic model instead yields QZ_ 12 U_ 1 CZ_ 12 , which is a 
two-qubit matchgate modulo local operations. For exam- 
ple, attempting = H_ x produces a maximally entan- 
gling gate that is locally unitary-equivalent to the match- 
gate Ghh [26] . In order to obtain a true local operation, 
required for universal MBQC when restricted to com- 
putational basis measurements, one must include some 
kind of fermion-fermion interaction. In the approach 
above, an additional hopping Hamiltonian with an on- 
site interaction yields an explicit CZ_ 12 which can cancel 
the effective CZ 12 arising solely from fermionic statistics. 
Because one can construct the SWAP gate from combi- 
nations of Hadamards and CZ gates, this approach is 
the measurement-based analogue of the fermionic circuit 
model found previously [221 ] . 

V. DISCUSSION AND CONCLUSIONS 

Through the 'leapfrog' model of fcrmions hopping 
across one another in a lattice, statistical antisymmetry 
and entanglement are shown to be equivalent quantities 
in fermionic systems. As a result, the ground state has 
maximal entanglement between all neighboring qubits, 
which is inherently insensitive to decoherence as long as 
the temperature remains lower than the energy gap. In 
principle, such a resource could be used for a variety of 



distributed quantum information protocols. Yet, univer- 
sal MBQC using this resource is impossible under a realis- 
tic model where only computational basis measurements 
are allowed. A universal set of gates can be obtained 
only when an additional local Hamiltonian is added that 
allows the fermions to explicitly interact through colli- 
sions. Note also that any double- well sites with zero or 
two fermions behave as defects (holes) in the cluster state; 
universal quantum computation is still possible as long as 
connectivity of the state exceeds the percolation thresh- 
old m. 

The results suggest that the ground states of inter- 
acting fermionic systems could generically be resources 
for quantum information processing within the MBQC 
framework. One requires only that the constituent 
fermionic wavefunctions have the opportunity to overlap 
one another, and interact in the process. Because these 
criteria arguably apply to almost all physical systems, 
discovering other suitable candidate Hamiltonians might 
be relatively straightforward, though finding convenient 
encodings for the quantum information might still pose 
a significant challenge. 
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